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Abstract 



A problem of constructing excited state swave functions of the discrete 
spectrum of completely integrable quantum systems is considered. Recur- 
rence relations defining wave functions up to the normalizing constant are 
obtained. 
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1. Introduction 

Explicit solutions of classical and quantum completely integrable sys- 
tems still attract much attention[l-5]. 

By now a lot of completely integrable systems with TV degrees of 
freedom have been discovered [3-5]. However for few of such systems ex- 
plicit solutions have been found. Some methods use additional constants of 
motion to reduce and simlify corresponding eigenvalue problems. 

In [6] we proposed a method of finding descrete spectrum energies and 
wave functions. This method doesn't use the fact of existence of additional 
constants of motion. In [7] we applied this method to the quantum systems 
with the following Hamilton function 

i=l ^ ' i>j 

with pair interaction potentials 

F(a;) = — L e = ,l 
sn z x 

External field for e = 1 has the form 

.... , MM" 1 ) 1 A(A + 1) 1 

^ = "2 ^ " "2 ^ A> ° (2<l) 
whereas for e = it takes the form 

W(x) = 2A 2 {e 4x -2e 2x ) (26) 

The first step of the method proposed in [6] is the transformation of vari- 
ables tj = (chxj)~ 2 , j = 1,...,N, which enables us to reduce the order of 
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sigularities in Schrodinger equation Hifj = Eifj. Then we perform furhter 
transformation to variables symmetrical in tj , j = 1 , . . . , N 

D L ~ L d 



ai = JJ^-, a z = -— — ai, D = ^ — , I = 1, AT 

j=i j'=i J 



which enables us to convert Schrodinger equation into the form convenient 
for solving it explicitly 

+ H 2 ((3)]k = [E-e(j3)]K 

N l-l 1 

= 2 ^2 [^2^2((l + m-2T)-is)a T a l+m _ T _„- (N -m + 2)a m _ 1 a l 

l,m=l t=1 v=0 

d 2 

-(N - I + 2)a m ai-! - (N + m + l)a z a m ]— — — (3) 

oaiddm 

N 

-^[a z _i(7V-/+2)(4/3+2^+3+2a(7V+/-3))+2a;(n-/+l)(2/3+l+a(7V+/-2)] 



1=1 



^ oaida 



l,m=l 
N q 

+ V a f — (4/3 + 2ju + 3 + 4a(7V - 1) + ) 

Here parameter (3 and e((3),p((3) are to be determined from normalizing 
condition. Wave functions are polinomials in symmetrical variables aj : 

N 

iP(x u x N ) = K (a u ...,a N ) H ^(1 - z]fl[(z^ - 



j=i j>fc 



= t/iXj and 



n 



K(a u ...,a N )=^2 c ju..,3Na J i-a J N ( 4 ) 

p=0 0<ji,...,jjv<JV 



Energy levels are determined by a set of TV integers ji , . . . , jjy , < j\ , . . . , Jn < 
n,ji + ... + jpj = p < n , where n - is the highest powerof the polinomial n. 
One can conclude from the normalizing condition that the upper boundary 
for n is 

n<^—^-a(N-l), e=l 

n< A- a(N - 1) - i e = 
Energy levels are the following 

TV 

-E h ,..., jN = -e(j3) - 4 J2 l Jdm + 2(iV + l)p(/3)(p(/3) + \ — n — 2n — aN) 

l>m 

N 

+2^2lj l [a{2N + l) + 2n + /i-X-j l -al} (5) 
z=i 

In [7] we only considered the simplest case n = 1 . Here we find wave 
functions of excited levels with n = 2. 

2. Poschl- Teller potential 
Expression (2) is a generalized Poschl- Teller potential. Having de- 
manded wave functions to have the following form 

N N 

«(ai, a N ) = ^2 °i a i + + X] bikdidk (6) 

Z=0 Z<fc 

and taking into account the normilizing condition we come to 

(3=^-a(N-l)-2 

p((3) = 1 - 2A (7) 
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m = -y [y^ 2 - 1) + (A - M - «(at - 1) - 4 ) 2 ] 



Substititing (5), (6) and (7) into (3), we come to a system of - 2 ' linear 
algebraic equations with zero determinant. Unlike the case n = 1, there is 
a contribution from the first term in Hi{(3) with second order derivatives 
in a,j. This leads to coefficients q, bik being dependent on all non-zero 
coefficients c m , m > /, 6 mn , m > Z,n > k. Performing induction one 
can obtain recurrence relations for coefficients in (6). As a result we come 
to 3 different sets of recurrent relations for different meanings of quantum 
numbers j 1 ,...J N . 

1. In this case ji = 26i s , s = 1, ...,N, p = 2; s' = 2s — N — 1, n' = 
2s — j, K = (N — s + 1)(A — fi — 2 — a(N — s)); 6 SS = <i s , where <i s is an 
arbitrary constant to be determined from the normilizing condition, and 



where Fjki are defined below. 

2. In this case ji = 5i s + Sin, s = 1, N — 1, p = 2; s' = s — 1, n' = 
AT, K = fi — A + 1; 6 s jv = <i s , where d s is an arbitrary constant to be 
determined from the normilizing condition, and 



bji = 0, j + l> 2s; Cj 



0, j > s 




s 



(8) 



bji = 0, j + I > s + AT; Cj - 



0, j > s 7 



Other relations have the form 



1 - 2A 
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b s -mN = — !— {-2[(JV - s + m + 1)(2A - 3 - 2a(N - s + m) + l]6 s _ m +i7v 



m 1 



+4(7V - s + m) ^2 ^2 b s - m+k N-k+r + (1 - 2A)c s _ m }, m = 1, s - 1 

fc=l r=0 

In this relation summation is performed under the condition 2k < N — s + 
m + t. 

It must be stressed here that Hamiltonian (1) is self-conjugated if and 
oly if (2A — 1) > 10 + 4a(7V — 1), i.e. the coefficients cj and bjk are not 
independent. 

3. In this case ji = 5i s + 8iN- m , s = 1, iV — 2, m = 1, iV — s + 
1 , b s N-m = d s N-mi s' = s — m — 1, n' = N + s — m — j, , where d s is an 
arbitrary constant, and 

bji = 0, j + I > s'; Cj = 0, j + l>N + s — m 

Coefficients 6 jn / are determined by (8). Suppose 

fljw = 2 Z n [(« - J + k - r)(A - n - 3 - a(27V — s — j + k — t)) — K], 

T=l 

fc-Z 

= 2k ~ i n t( s - - ^ - 3 - a t 2N - s - j + !)) 

T = 

+ (7V - n' + fc - r + 1)(A - fi - 1 - a(7V - n' + k - t)) - K], 

S j = R jll + N-j + l, 
Qj = Rju + 2(7V - s + 1)(A - - 3 - a(7V - 3)) + 2(2A - 2fi - 7). 

Then the common part of the three recurrence relations sets is the following: 

(N — s' + k)\ 
Cs '~ k+l = r ((2A - 5) /(2a) -N + s'-k) [{ ~ 2a)k 
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r((2A-5)/(2a) -N + s') 

(N-s')\R a , kk Cs ' +1 

A r((2A-5)/(2a)-iV + s'-fc + Z-l) 

+ ^ } (N-s> + k-l)\R s , kl 

s-s'+k-r-l 1 

^ ^ 6s^-fc+;+ -+rAT- -+l(l — &rl<fyfc]; & = l,...,s' 
cr=l T =0 

In this relation summation is performed under the condition 2<j < N + 1 + 
k — s f — I — T. 

k (TV-n' + fc + l)! r r((A-2)/a-jV + nQ 
Djn/ - fc_l aj r((A-2)/a-7V + n '_/^ (7V- n ' + l)!F jfcl °' n ' 

-2^-,.,) v y^c-^)* r ((A ~ 2)/Qg - ^ + ^ - o r 

(2(N — j + 1)(A — 2 — a(N — 1)) +n f — I + j)b j+ln ,-i 

s 1 

-2(n-l-j) ^2 ^2b an >-i- a+T+j (l - 5 Tl 8 aj+1 )}}, k = l,...,n'- j + 1 

a=j+l t=0 

In this relation summation is performed under the condition 2a < n' — I + 

1+T + j. 

_ 2(N-j + l)(2\-3-2a(N-j)) u 
d jj Ojj+i- 

3. Morse potential 

One can come to a a case of molecular Morse potential by a limiting 
procedure Xj = q, q — > oo. As a result potential (2) turns into (2b), and 



A - ii -> 2 A + 1 
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(9) 



However it's easier to find wavefunctions directly from (3) and perform the 
calculations similar to those in section 2. Substituting (9) into (5) and 
(7), one can derive energy levels in Morse potential. Taking into account 
p{(3) = — 8A recurrence relations for the coefficients in (6) turn out to be 
the following 

(— 2A) k 

c s ,. k+1 = (N - s> + kn j^^- c, +1 

k ( — 2A) l ~ 1 S-s'+k-T-l 

+ 4 E (jV _^ + fc _ ;)!i } s , M E b s ,. k+l+<r+rN . a+1 ], k = i,...y 

In this relation summation is performed under the condition 2a < N + 1 + 
k-s'-l. 

^, k = (^-»' + t+ l)! {r JS_ 

-2E (y _ir;f|\ )!fj j (^-i + i)w-, 

s 

-2(n'-l-j) ^ &an'-Z-a+j]}, = 1, n' - j + 1 
o-=J+l 

In this relation summation is performed under the condition 2a < n' — I + 

& S -miV = 77^— {[(# " 8 + m + 1)(1 " 4A) - I]6 s _ m+ 1TV 

m 

+4(7V - s + m) ^2 b s - m +kN-k ~ 4:Ac s - m }, m = 1, s - 1 
fe=i 
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In this relation summation is performed under the condition 2k < N — s+m. 

_ ±A{N-j + \) 

°33 - + 

It is sufficiently to perform A — fj, — > 2 A + 1 in the other relations. 

Therefore we have found recurrence relations which determine coffi- 
cients cj and bjk for arbitrary N for Poschl- Teller and Morse potentials. 
Wave functions are determined by means of (6). 

Let us stress here that at < A — \x — 2a(N — 1) < 6 and 1 < 2A — 
2a(N — 1) < 7 wave functions found together with wave functions for n = 1, 
which were found in [7] , give the whole discrete spectrum of the systems in 
question. 

As an example we have found four particle wave function of one of the 
excited states in Morse potential (see Appendix). 

The results obtained could be used for testing of various approximate 
methods of solving many particle quantum system as well as in supersym- 
metrical generalizations of such systems. They could also be useful for 
investigations of systems close to integrable. 

Let us note in conclusion that wave functions of higher levels with such 
values of parameters A and a that n > 3 can be found in the same way. 

Appendix. 

Here we write down explicit form of the coefficients cj and bjk for the 
14th exited state for four-particle system in Morse potential. Parameters A 
and a satisfy inequality 1 < 2A — 6a < 7. Discrete spectrum consists 
of 15 levels: the ground state, 4 levels with n = 1 (p = 1) and 10 levels 
with n = 2{p = 2). The energy of the level in question corresponds to 
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ji = 8 13 + 5/4 and is equal to 

E = -2(4A 2 - 12aA - 16A + 14a 2 + 30a + 19) 

Suppose 

F(A, a) = (1 - 2A)(12A - 5a - 23) (4A - 3a - 6) (4,4 - 3a - 5) 

T(A, a) = 2(24A 3 - 104A 2 + 42aA + 116,4 - 3)(4A - 3a - 5) /F(A, a) 

X(A,a) = 8A(4608A 4 - 5376aA 3 - 15720A 3 + 1035a 2 A 2 

+9504aA 2 + 13144A 2 - 42aA - 116A + 3) /F(A, a) 

y(A, a) = 3(73216A 5 - 248565A 4 - 85632aA 4 + 150328aA 3 

+149893a 2 A 3 + 204520A 3 + 1888aA 2 + 2285A 2 

-624A - 126a 2 A - 606aA + 9a + 15)/[(8A - 6a - 9)F(A, a)] 

Z(A,a) = 64A 2 (32A 2 - 24aA - 36a - S)(UA - 10a - 9)/(8A - 6a - 9) 

296A - 366a - 421 



N(A,a) = 
M(A, a) = 



8,4 - 6a - 9 
2A 2 [3Z(A, a) - N(A, a)Y(A, a)] 



(5A-6a-7)(3A-4(a + l)) 
Then the coefficients of (6) depend on the potential depth A and pair in- 
teraction constant a in the following way 

^ = 2 { 3A-Aa-^ a) + 2^1 ) - (2A-2a-mA-l-6r ) d 

6 

C2 "4A-3a-6 d 
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4A 

bn=M {A a) ZA _ i{a + l) d 
bi2 = M(A,a)d 

-4A 

bis = TT—a v\Y( A > a ) + Z(A,a)]d ■ b 14 = Y(A,a)d 

oA — ba — 7 

12A 

622 = 8A -6a-9 X ( A ' ; 623 = 

8A 

624 = T(A, a)d ; 6 33 = 2A l d 5 hi = d 
Here d is an arbitrary constant to be fixed by normilizing condition. 
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